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The energy and structure of dilute gases of hard spheres in three dimensions 
is discussed, together with some aspects of the corresponding 2D systems. A 
variational approach in the framework of the Hypernetted Chain Equations 
(HNC) is used starting from a Jastrow wavefunction that is optimized to 
produce the best two-body correlation factor with the appropriate long range. 
Relevant quantities describing static properties of the system are studied as 
a function of the gas parameter x = pa d where p, a and d are the density, 
s-wave scattering length of the potential and dimensionality of the space, 
respectively. The occurrence of a maximum in the radial distribution func- 
tion and in the momentum distribution is a natural effect of the correlations 
when x increases. Some aspects of the asymptotic behavior of the functions 
characterizing the structure of the systems are also investigated. 

PACS numbers: 03.75. Hh, 05.30. Jp, 67.40.Db. 



1. INTRODUCTION 

The study of dilute systems has become a subject of major interest since 
the achievement of Bose-Einstein condensates in low-density atomic gases 
confined in harmonic traps. The gas parameter x = pa d where p, a and d are 
the density, s-wave scattering length and space dimensionality (2 or 3 in the 
cases under study), govern the behavior of the gases in the dilute regime. In 
the 3D case, an universal dependence of the energy on the scattering length 
alone has been proved to hold up to values x ~ 1CP 3 , while deviations appear 
at higher x. 

In three dimensions, low density expansions can be derived in the frame- 
work of standard perturbation theories. Infinite sums of ladder diagrams can 
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be carried out and result in the well-known Lee and Yang expansion of the 
energy per particle^ 
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measured in units of Ti 2 /2ma 2 . 

In two dimensions, logarithmic divergences in the scattering length in- 
troduce additional difficulties in the derivation of similar expressions. Still, 
Schick and Lieb 12121 proved that at low x, the energy per particle of a 2D gas 
of hard-core bosons satisfies the inequality 
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in the same units 

In this work an analysis of the energy and other properties of 2D and 
3D gases of Hard-Sphere bosons is presented ™. The starting point is a 
variational wavefunction of the Jastrow type 



*o(ri,r 2 , ...,r n ) 



n fM 

l<i<j<N 



(3) 



which is known to accurately account for most of the energy and structure of 
homogeneous Bose systems at low densities. The minimization of the energy 
corresponding to the Hard-Spheres potential 



V(r) 




(4) 



produces a suitable wavefunction that can be used to calculate the energy as 
well as other, relevant quantities describing the ground state of the system. 



2. HNC FORMALISM 

The energy per particle of an homogeneous gas of Hard-Sphere bosons 
described through a Jastrow wavefunction can be expressed in terms of the 
two-body correlation factor f(r) 

E = -\p J drg(r)^-V 2 f{r) (5) 

and the radial distribution function g(r), 

N(N -I) J cfr 3 dr 4 • • • dr N | * | 2 
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One can formally express f(r) as a functional of g(r) and solve the 
Euler-Lagrange problem 5E[g]/5g(r) = 0. This leads to a set of equations 
which are easily expressed in terms of the static structure function, S(k) = 
l + p/dre ik ' r (s(r) - 1). These can be solved and yield 

S(k) = t{k) , (6) 

y/t*(k)+2t(k)V ph (k) 

with t(k) = h 2 k 2 /2m and V p h(k) the Particle-Hole interaction P. In the 
HNC/0 scheme used throughout this work, elementary diagrams are dis- 
carded due to their small contribution at low densities. Then, for a system 
of Hard Spheres and in r space 

h 2 ^ r-r |2 



V P h(r) = - | VJg(r) \ 2 + [g(r) - 1] u//(r) , (7) 



where the fe-space induced interaction coj(k) reads 

(h , 1 [25(fc) + l][5(fc)-l] 

Once with S{k) or g{r) one can get f(r) and compute the momentum 
distribution n(k) 

n(k) = (2n) d pn S(k) + J dre^[ Pl (r) - pn ] , (9) 

where no is the condensate fraction and pi(r) is the one-body density matrix 
/ s AT JdT2---dT N ^ Q (Ti,r2,...,T N )^o(ri',T2,...,r N ) , 

Pi(m>) = N — — — r-— j2 • io 

J dridr 2 • • • dr N | V \ z 

The momentum distribution is normalized as 1 = J dkn(fc)/(27r) <i p and 
its second /c-weighted moment provides the kinetic energy per particle. 

Notice that for Hard Spheres the potential energy is zero because the 
product g{r)V{r) vanishes, and therefore the kinetic energy coincides with 
the total energy of the system. 



3. RESULTS 

The energy per particle of the gas of Hard Spheres gas is plotted as a 
function of the gas parameter x in the left panel of Fig. Q where the Euler- 
Lagrange solution is compared with Diffusion Monte Carlo results, the low 
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Fig. 1. Energy per particle of the gas of Hard Spheres in 3D and 2D (left 
and right panels). Solid circles, stars and open squares in the left stand 
for Euler-Lagrange, HNC/0 and Diffusion Monte Carlo results (the latter 
taken from Ref. 7. The dashed line represents the low-density expansion of 
Eq. (^Q). Solid circles and open squares in the right are Euler-Lagrange and 
HNC/0 results from the short-ranged/s/j(r) of Eq. (|12j) . 



density expansion of Eq. (^Q) and the HNC prediction obtained from a short- 
ranged two-body correlation function fsR.(r). This last function minimizes 
the lowest order in the cluster expansion of the energy of the homogeneous 
gas of Hard Spheres with a healing condition at a distance d, taken as a 
variational parameter^ 

( r < a 

fsR(r) = < d sin[Jf(r-q)] , r , . n , (11) 
{ r sm[K(d-a)} r > a 

where K fulfills the equation cot[K(d — a)] = (Kd) . As can be seen, 
differences between the different approximations are not significant up to 
x ~ 0.001. The leading term in the expansion, E(x) = -iirx, is not enough to 
correctly account for the energy of the system even at lower values of x, while 
the inclusion of the second term, which introduces additional dependence on 
x, improves the agreement. The influence of the optimization on the energy 
is rather small and the energy is dominated by the short range structure of 
the potential. 

The scaled energy per particle of the gas of Hard Spheres in 2D is shown 
in right panel of Fig.^ Euler-Lagrange results are plotted in solid circles and 
compared with the energy produced by a short range two-body correlation 
with a healing distance, as in the 3D case. Geometry constrains change the 
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Fig. 2. g(r) for the 3D (left) and 2D (right) gases of Hard Spheres as a 
function of the gas parameter. Solid line: x = 1CT 1 , dot-dashed line: x = 
10~ 2 , long-dashed line: x = 10 -3 , short-dashed line: x = 10~ 4 . 



shape of fsR,(r) which now reads 

[ r < a 

f{ r ) = \ Y (\a)MXr)-Jo(\a)Y (\r) r . n (12) 
L Yo(Xa).Jo(\d)-J (\a)Yo(Xd) ' ^ " 

where A is a Lagrange multiplier satisfying the condition Ji(Xd)Yo(Xa) = 
Jo (Ad) Yi (Act). Energies have been divided by the leading term in Eq. @ to 
check convergence towards the limiting case. The agreement of the energy 
with the corresponding low-density expansion is similar to the 3D case when 
only the first term E{x) = Airx is considered, while deviations start to be 
significant at x ~ 0.001. In any case, the singular behavior of the energy 
seems to be entirely taken into account by the logarithmic term, and thus the 
remaining deviations are expected to be successfully modeled by a smooth 
function of x that approaches zero when x decreases. 

The radial distribution functions of the 3D and 2D gases of Hard Spheres 
are shown in Fig. |2 for several values of x (left and right panels, respectively). 
In both cases g{r) develops a peak at x = 0.1 that can not be resolved at 
lower values of a;, a clear signature of the effect of correlations. At low 
x's, g(r) is a monotonically increasing function of the distance, approaching 
faster and faster the asymptotic limit g(r) — > 1 with the density. This is 
more notorious in the 2D case. In fact, this is not surprising if one notices 
that in 3D and at large distances g(r) ~ 1 — a/r 4 , while in 2D this behavior 
changes to g(r) ~ 1 — 6/r 3 , with a, b constants depending on the velocity of 
sound in the medium. 

The final quantity analyzed is the momentum distribution. The product 
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Fig. 3. Momentum distribution of Hard Spheres in 3D and 2D (left and right, 
respectively). Solid and dashed lines in the left, x = 0.05 and x = 0.08. Solid 
line in the right, x = 0.01. 

kn{k) is shown for the 3D and the 2D cases in Fig. |31 At low x, the mo- 
mentum distribution is a monotonically decreasing function of k, a behavior 
already shown by the Bogoliubov approximation at any density. When x 
increases, kn(k) develops a peak at low k, which can be taken as a gen- 
uine effect of long-range correlations. This peak becomes more pronounced 
and shifts to the right with increasing x. While this maximum is clearly 
visible already at x = 0.01 in 2D, it only shows up at x > 0.05 in 3D. As 
far as the condensate fraction uq is concerned, the Bogoliubov estimation 
in 3D, Uq = 1 — (8/3)1/^7^) predicts values that are slightly higher than 
the ones obtained from the solution of the Euler-Lagrange problem, yielding 
Uq = 0.850 against no = 0.801 at x = 0.01. Larger discrepancies arise at 
higher values of x. Similar results hold when the EL condensate fraction is 
compared with Schick's prediction, Uq = 1 + l/ln(x). 
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